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Abstract 

We calculate the QCD corrected effective Hamiltonian for BB-Mixing 
in heavy quark effective theory including corrections of the order 
^Qcn/'^b- The matrix elements of the subleading operators are es- 
timated using the vacuum insertion assumption. We show that the 
major part of the subleading corrections may be absorbed into the 
heavy meson decay constant Jb] the remaining corrections are only 
due to QCD effects and give an enhancement of AM of 5%. 
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1 Introduction 



Oscillations between particle and antiparticle have been observed first in the 
neutral Kaon system and later also in the neutral B meson system P, |]. 
Such oscillations are predicted by the standard model and proceed through 
the so called box diagrams depicted in fig. 1. Since the top quark as well as 
the W boson are heavy compared to the B meson it is convenient to describe 
the oscillations in an effective theory, where the top quark and the W boson 
are integrated out. 

We shall concentrate in this paper on the oscillations in the neutral B 
meson system. The effective Hamiltonian, which is obtained after integrating 
out the top quark and the W boson, consists of only one operator 



= -I\V,lVa\'mii> ( I 0„ (1) 

where 



Oo = (dLl.b) (dLYb) (2) 

{(11 denotes the left-handed component of the d quark field) and $ is a 
function which arises from integrating out the top quark and the W boson 
at the same scale 0] 

1 9 3 3 x^lnx , . 

^"^^ " 4 ^ 4(1 -a;) " 2(1 - s)2 " 2 {1 - ^ ^ 

To evaluate the mass shift AM relevant for the oscillations, one needs to 
calculate the matrix element of this effective interaction between a neutral 
B meson and its antiparticle; it is given by 

AM = ^(fi°|/7eff|5^) 



F 



27r2 



(1;) ^m^} (4) 



where fs is the B meson decay constant defined by 

{0^,^,h\^{p)) = tfBP,. (5) 
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The coefficient ?7qcd contains the short distance contribution of the QCD 
corrections which is calculable in perturbation theory. It depends on the fac- 
torization scale fi at which the short distance contributions become separated 
from the long distance ones; the non-perturbative long distance contributions 
are contained in the "bag parameter" Bb, which depends in such a way on 
/i that the matrix element of the Hamiltonian is scale independent . 

As an estimate for the order of magnitude of the effect the so called 
"vacuum insertion" assumption has been frequently used in the past, which 
corresponds to the case Bb = 1- However, since Bb is scale dependent this 
assumption has to be supplemented by a statement about the scale where it 
is assumed to be valid. 

The short distance contribution from scaling down to a scale < fi < 
Myy has been calculated in leading logarithmic approximation some time ago 
[H] and yields 



-6/23 



^QCd(/^) = f^#T ) mb<fi<Mw, (6) 



which gives a correction factor of r^qco ~ 0.85 at the scale of the b quark. 
The corrections induced by subleading logarithms wich are to be considered 
in conjunction with the 0{as) matching conditions at the W scale have also 
been calculated p. 

If one could get some information on the parameter i?^ at a scale of the 
order of the B quark mass, one could in fact predict the amplitude of the 
oscillations as a function of the CKM matrix elements and the top quark 
mass. However, at the present time the only non-perturbative information 
available is from lattice methods 0, which gives some information on the 
parameter JbBb at low scales^. Thus it is desirable to scale further down 
to these very low scales. This is possible by switching at the point fi = rrib 
to heavy quark effective theory |p, in which the b quark is integrated out 
and replaced by a static color source. The leading logarithmic result for the 
scaling /i < m^, i. e. the resummation of logarithms (a^ ln(mf,//i))", has been 
calculated in The result is 

VqcdW = \~^^ I VQCDinib) mc<n<mb (7) 



^In fact, the lattice results indicate that the "bag parameter" Bb at these low scales 
is indeed close to the value Bb = 1 used in the vacuum insertion assumption 
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/ , X X 12/27 

^QCD(yU) = T ^QCd("^c) ^<me. (8) 

This result is indeed very peculiar, since the anomalous dimension of the 
four fermion operator Oq is the sum of the anomalous dimensions of the 
two left handed currents involved. In fact, if factorization like e. g. vacuum 
insertion assumption takes place, one would expect this behaviour of the 
anomalous dimension to all orders. However, the QCD scaling below rrib 
has been calculated even to subleading order |jlO| and it seems that this 
factorization does not hold for the subleading result. 

Heavy quark effective theory is a systematic expansion in inverse powers 
of the heavy quark mass. The purpose of the present paper is to extend the 
existing calculations of the scaling below by including terms of the order 
l/ruf, into the effective Hamiltonian. The motivation is twofold. Firstly, 
there are indications from the lattice that the l/rrih corrections could be 
large, even for the b quark. Secondly, this is the first attempt to deal with 
the l/rrib corrections in a non-leptonic process, of which the Hamiltonian for 
BB mixing is certainly the easiest to deal with. 

In section 2 we shall discuss an appropriate operator basis for the l/rub 
corrections to BB mixing and calculate the anomalous dimension matrix and 
the Wilson coefficients at some low scale fi. In section 3 we shall estimate 
the 1/mfe effects by using the vacuum insertion assumption and the results 
on the 1/m corrections for the heavy meson decay constants obtained from 



QCD sum rules [|TT], |T| . 



2 The Operator Basis in the Order I /nib 

To leading and subleading order in 1/m, the effective Lagrangian for a heavy 
quark moving with velocity v is given by (we use = — igA'^^T'^) 

L = htiv-Dht + ptK + Kpt (9) 

where the transverse derivative D — -^v-D is denoted as D-^, the heavy quark 
field is defined by 

= e'"^''-"^4^6(a;), (10) 



3 



and the sources and for the heavy quark and antiquark fields with 
velocity v have been retained. Since there is also a heavy antiquark involved, 
one has to use also the effective Lagrangian for an antiquark which can be 
obtained from (pi) just by replacing v with —v and /i+ with 



1 - 



-h{x) 



(11) 



along with the proper replacements of the sources: ^ , . 

Using the quark-type parameterization for one factor in the effective 
Hamiltonian and the antiquark-type parameterization in the other fac- 
tor (using the nomenclature of [jl3[) one obtains the effective Hamiltonian 
density at the matching scale p = mb in terms of heavy quark fields 



eS 



mi 



(12) 



X 



nib 



ai{mb) X^{mb) + ^Ci(mb) C+(mfc) 



i=l 



At the matching scale there are two local operators 



(13) 



and three nonlocal operators , which are time-ordered products originat- 
ing from the 1/m piece of the Lagrangian: 



Of 

ot 



d^xT[O'Q{0) CK{^v ■ Dfh+ + h-{iv ■ Dfh-){x 
d^xT O'oiO) CK^cr^'G.X + h-\^''"G,,h-){x) 



(14) 



The coefficients of the l/mf, terms are 



ai{mb) = |?7qcd("^6), 
C2(mfe) = Csinib) = -ci{mb) = ^r]QCD{mb), 



(15) 
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In order to calculate the leading logarithmic QCD corrections to this 
effective Hamiltonian, one first has to set up a basis of operators bearing the 
correct dimension and quantum numbers. Since there are three independent 
momenta in the transition, and five independent gamma matrix structures 
possible, one arrives at a basis of 15 local operators: 

Pi = {dLiv ■ Dh^) {dLK) 

Ps = {dLiDX){dLrK) 
Pa = {dLi^v-Dh^){dLl^h-) 

P5 = iex^^upv\dLiD''rK){dLYK) 



Qi = (dLh+) (dLW ■ DhZ) 

Q2 = {dLlX){dLiD^h-) 

Qs = (dLK) (dJI^h-) 

Qi = {dLli,K){dLl^w ■ Dh-) 

Q5 = iex^.y{dLrK){dLiD''YK) 



(16) 



Ri 
R2 
R3 
R4 

R5 



iv ■ d{dLhl) {(IlK) 
id^,{dL^u,h:^)] (diK) 
idf^idiK)] (diYK) 



id^{dLYh+)] {dLYK) 



Operators containing 75 or a matrices can be eliminated by the projection 
operators implicit in the d^ and fields. 

These operators can mix with each other through the diagrams of fig.2. 
In addition, there are also the diagrams of fig. 3 which introduce mixing of 
the nonlocal operators Of with the local operators. Note that there is no 
mixing in the other direction since the local operators do not require the 
nonlocal ones as counterterms. 

One can employ the symmetries of the effective Hamiltonian to simplify 
the mixing matrix. Firstly, the lowest order operator as well as the 1/mb cor- 
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rections in (|T2]) have the property of being symmetric under Fierz transfor- 
mations, i. e. after exchanging the two hght quarks one recovers the original 
operator by rearranging the indices of the gamma matricesQ Therefore it is 
useful to switch to an operator basis {Xj, Yj} of definite Fierz parity: 

J^Xi = Xi and TYj = -Yj. (17) 

In the case at hand, there are seven linear combinations of P, Q, and R with 
positive and eight combinations with negative Fierz parity. 

Secondly, the operators can be divided into a class of operators {X^^, Y^^} 
which are symmetric under the exchange v —v, and another class of 
antisymmetric operators {X~ ,Y~}. It turns out that there is no mixing 
between these two sets. Since the operators introduced by tree-level matching 
(|T^ have positive parity both under Fierz and v —v transformations, the 
X~ and Y^ operators can be neglected altogether. 

Thirdly, one can eliminate all operators which vanish because of the equa- 
tions of motion 

v-Dh^ = 0, I^dL = (18) 

since they do not contribute to physical matrix elements and do not mix with 
operators that do contribute. 

The anomalous dimensions of the nonlocal operators is obtained simply by 
adding the anomalous dimension of the lowest order effective Hamiltonian Oq 
and the anomalous dimensions of the 1/m terms in the effective Lagrangian 
[Q. Since the former is just twice the anomalous dimension of the involved 
current, this pattern of factorization is found again for the nonlocal part of 
the recoil corrections to the four fermion matrix element. 

Evaluating the one loop diagrams depicted in fig. 2 and 3 we obtain 
the result summarized in tab. (tab. 0): The mixing matrix of the relevant 
operators separates into two blocks. The chromomagnetic moment operator 
mixes with the Fierz symmetric operators 

X+ = P2 + P3-P5 + Q2 + Q3-Q5 

X+ = ARi + Ri (19) 

the first one being the operator that already occurred in the matching condi- 
tions (P^. The kinetic energy operator O2 requires as counterterm a Fierz 

^ For the Ri operators which arise in the mixing, in order to obtain the Fierz transform 
a partial integration is also necessary. 
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Table 1: The anomalous dimension matrix of the relevant operators. A factor 
as/12'K has been omitted in all elements. 



antisymmetric operator which is renormalized multiplicatively 



(20) 



With these definitions, the effective Hamiltonian at scales < mb is given 
including the 1/mb terms by 



eff 



mi 



r/QCD(/^)Oo(/i)+ —[c2{^^)Ot{^^){2l) 



+ cM Otif^) + ai(/i) X+ifi) + aM X+ifi) + 6i(/i) Y+{^^)}' . 
The array of Wilson coefficients {ck,ai,bj} is the solution of the renormal- 



ization group equation 




H 7 

127r ' 



( 

a 







(22) 



where 7 is the anomalous dimension matrix as listed in tab. [1|. 

The entries in the anomalous dimension matrix have been calculated in a 
general covariant gauge and have been found to be gauge independent, as it 
should be the case. However, there is a gauge dependence in the counterterms 
of the Oi operators corresponding to operators which vanish according to the 
equations of motion ([l8|) , indicating that these terms are unphysical. 

The solution of fE2|) can be written in the form 




exp 



7Tln 



2(33 — 2nf) as{mb) 



( c{mb) 
a{mb) 
-I \ b{mb) 



(23) 
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where rif = A for /i > mc, and n/ = 3 for /i < nic- With the initial conditions 
at II = nib taken from (|15|) the solution (^31) reads 

C2(/i) = C(/^)'' 
C3(/i) = C(/U)^ 



with 



= -IcM'- 4(413 + VM9^ .,_^,, 
Uf-; 649(-27 + ^649) 

+ 4(413 -v^) ,«,^V, 
649(27 + ^649) 
a.(,) ^ 11^(^)3 _ 40(295 + llvg9) ,,_^,, 
^'^^ ir^^^ 649(27+^649) 

^ 40(295 - llv/649) (4,^^^/4 
649(-27+ V649) 
&i(-") = -C(^)' + C(/i)'° 



\as{mb) J \as[mc) ) 



for /i < TTT-c. Inserting the values 

mft = 5 GeV, = 1.5 GeV, Aqcd,4 = 0.2 GeV (26) 
and the one loop approximation for we obtain at /i = 1 GeV 

(c2,C3,ai,a2,6i) = (0.60, 0.47, 0.52, -0.16, -0.15) (27) 
together with the value 77qcd(1 GeV) = 1.20. 



3 Estimate of the Matrix Elements 

In order to state numerical results, one has to calculate the matrix elements 
of the operators as given above, evaluated between eigenstates of the lowest- 
order HQET Langrangian since the mass dependence resides now in the 
coefficients. As long as a non-perturbative calculation is not available, one 
has to rely on model assumptions. We shall estimate all matrix elements 
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using the vacuum insertion assumption. As pointed out in the introduction 
we expect this assumption to at least give the correct order of magnitude 
since the parameter Bb calculated from the lattice is indeed of order unity. 

Using the vacuum insertion assumption for the matrix elements of the 
effective Hamiltonian we shall encounter the same matrix elements as in the 
discussion of the heavy meson decay constant Jb which has been analyzed 
to order I/itlb using QCD sum rules |11, 12|. The meson decay constant 



fs and its leading recoil corrections may be parameterized in terms of four 
constants, which have been estimated in Both results for J'b [0! HI 



are 



consistent and we shall use the values given in [|TT| for the four parameters. 

The parameters describing the heavy meson decay constant are defined 
by the following matrix elements 



{0\dTh+\B\v)) 



\^^,idTht)\B\v)) 



lF(/.)Tr{rM+(t;)} 
^F{i2)Av^Tt{tM+{v)} 



I I d^x {0\T [idTh+){0) {h^{tDfh^){x)\ \b\v)) 
F{fi) Gi{i2)Tt{tM+{v)} 



d^x (0|T 



(rfT/.+)(0)|(/i+a^.G"-/i+)(x) 



\B"{v)) 



(28) 
(29) 

(30) 
(31) 



6F(/i) G'2(/i)Tr{rM+(t;)} 



The constant A is the difference between the heavy quark and the heavy 
meson mass; the other three constants F{fi), Gi(/i) and 6*2 (/x) depend on the 



renormalization point fi. In [|rT| these constants have been estimated at the 
low scale n ~ 2A. The values obtained from QCD sum rules for the four 



parameters are [11 



A = 500 MeV, F(2A) = 0.36 GeV^/^ (33) 

G'i(2A) = -0.5 GeV, G'2(2A) = -55 MeV 

where the effective value for Gi as discussed in [|ll| has been used. Fur- 
thermore, the matrix M+(f) appearing in eqs.( |2^ - |3T| ) is the representation 
matrix for a heavy pseudoscalar meson containing a heavy quark 



M+{v) 



-v^mi(l + 1^)75. 



(33) 
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Vacuum insertion for the lowest order operator O'q corresponds to the 
replacement 



4 



{B^\dL^^h-\0){0\dLrK\B 



(34) 



where the factor 4/3 is a colour factor: The light quark operators may be 
contracted in two ways with the light quarks in the mesons which corresponds 
to taking the Fierz transform. One of the possibihties is colour singlet and 
contributes with a factor 1, the other one is a combination of colour singlet 
and octet, of which only the singlet contributes. This yields in total a factor 
of 4/3 for a Fierz symmetric operator like O'q. Similarly, one would obtain a 
factor of 2/3 for a Fierz antisymmetric operator. 

To evaluate this the matrix elements involving mesons containing a heavy 
anti-quark arc needed. They are given by the same constants A, F{ii), 
and G2(/u) due to charge conjugation and are evaluated by replacing the 
matrix M~^{v) by 



I 



(35) 



Thus we obtain for the matrix element of O'r, 



{B'n\B")^-F'{t,)mn. 



(36) 



After vacuum insertion all the matrix elements of the local operators in 
the order I/uib Pi, Qi and Ri may be expressed in terms of -F(/x) and A. The 
non- vanishing matrix elements of the local operators are 



{B' 
(5° 



P2 



b') 



P,\b') = {B'\Q,\B 



R2\B . 



^0, 



12 



AmsF^ilj,) 



Rim - 



--ArriBF'ii,) 



-{B^'lRslB 
1 



--AmsF\^) 



6 



AmBF\n) 



(37) 
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where the two possible contractions of the hght quark in the four fermion op- 
erators have been taken into account. For the local operators with a definite 
Fierz parity, the result is 

2 _ 

— AmBF^(/i) 



(5°|X+|F^) 

(S°|X2+|#) 



(38) 



We shall now define the vacuum insertion for the non-local operators Oi~. 
For the piece of O2 containing the quark operators we define vacuum 
insertion by 



rf^x(50|T Oo(0)(/i+(zD)2/i+)(x) \b\v)) 



(39) 



where the colour factor of 4/3 is the same as in the lowest order piece Oq. 
Thus this matrix element may be expressed in terms of F{n) and Gi{n). 
Analogously we define vacuum insertion for the matrix element of the non- 
local operator and obtain 



(40) 



{B'\Ot\B') = 8mBF'ifi)G2ifi) 



while the matrix elements of Oi vanish due to the equations of motion of the 
heavy quark. 

Within the framework of the vacuum insertion assumption, we now can 
insert the expressions for the matrix elements (^Bp and (^) into the effective 
Hamiltonian (pTI): 



AM 



M. 

A 



VQCBif^) + — i4c2(/i) 



rrib 



A 



+ 24c3(/i) 



A 



- 2ai(/i) - 3a2(/i) + 61 (/x) 



(41) 
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With the numerical values taken from ( p7D and (0), appropriate for the scale 
11= 1 GeV, we obtain 



g|V,W.M*(^)F=(lGeV)(l.20-4.4A) 



' 'V^;,V^*.r-?$ ( 47r 1 0.51 GeV^ (42) 



Thus, with the parameters taken from QCD sum rule analysis, the 1/mb 
terms amount to a large correction of about -40% in this case. The most 
important l/m^ corrections arise from the nonlocal operators; furthermore, 
the fact that Gi is not well known |jTl| introduces a large uncertainty in the 
final result. 

However, the matrix elements of the nonlocal operators factorize in lead- 
ing log order in the same way as the leading operator O'q does, because the 
diagrams of fig. 3 which could lead to non-factorizable contributions only 
introduce local operators as counterterms, and the local operators in turn 
do not mix with the nonlocal ones. Therefore one can easily absorb these 
contributions into the pseudoscalar decay constant /b, if the contributions 
of subleading local operators to are taken into account properly. 

Explicitly, the square of the pseudoscalar decay constant is given in lead- 
ing log approximation up to order 1/m;, by |]TT], |K| 

Jl = :^C"(M) f 1 + 2^ + UC-'M^ - A |i + d(,)]) , (43) 

TUB V mb nib J 

where for /i < 

1/25 / / \ \ 1/27 



C(,)=(^) 1^1 . (44) 

\as{mb) J \as{m^) J 

The radiative corrections in the local term 

1 fi 

d{fi) = ^{C-\fi)-l-\nC\fi)) (45) 

remain very small for reasonable values of /i. In particular, d{l GeV) = 0.05. 

Comparing the coefficients of the lowest order operator and the nonlocal 
l/rUb terms in (|T) and (^) shows that these terms factorize and we have 

VQCBif^) = C2(/i) = C^^(^)^QCD("^b) 

C3(/i) = C^(/^)^qcd("^6)- (46) 
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This may be used to rewrite ( ^T]) by absorbing the large nonlocal l/nib con- 
tributions into 



M2 




1 + d{^i) (47) 

H 7 — -(-2ai(/i) - 3a2(Ai) - asl/w) - + &2(/u)) 

r]QCD{mb) 

This expression has the form of eq. (^) from which we may read off the bag 
factor 

BBimb) = 1 + 0.45— = 1.05 (48) 

rrib 

where /z = 1 GeV has been taken as the scale where vacuum insertion is 
assumed to be valid for the subleading matrix elements. 

In fact, the leading log result only implies that in the static limit Bsifi) 
is scale independent for A < < rub, where A is the scale where perturbation 
theory breaks down. However, the static value i?^'^* of -B^ is not fixed and 
we have simply set its value to unity. For an arbitrary value of iJ^*** our 
results for AM would be multiplied by this factor. In particular, eq. (^) 
also is multiplied by i?^*^* and thus we have calculated the l/rub contributions 
to -Bs as defined in (^. We note that these corrections can be calculated 
perturbatively. 



4 Discussion of the results 

We presented a complete calculation of the leading log QCD corrections 
to order l/m^ for the effective Hamiltonian relevant for BB mixing. The 
matrix elements have been estimated using the vacuum insertion assumption 
which reduces matrix elements of four fermion operators to a product of 
current matrix elements between the meson and the vacuum. The latter 
have recently been estimated using the QCD sum rule approach |11]. As far 



as Jb is concerned the l/rub corrections according to are indeed large; 



however, in our analysis we have absorbed these large corrections into fs- 
The remaining corrections are calculable subleading corrections to the bag 



13 



parameter which are of the order (as/'7r)(A/m;,) \n{mh/^) with /i = 2A. They 
arise only due to QCD effects and hence are smalL We find an enhancement 
factor for Bb of 1.05. 



We may also compare the results to experiment [IC]. Using a value for x 



X = AM ■ Tfe = 0.67 ± 0. 10 (49) 

we may study the possible values of |V^d| as a function of the top quark mass. 
Keeping in mind the above assumptions and using a value of fs obtained 
from the lattice 

fs = 205 MeV (50) 

we find 

0.008 < \Vtd\ < 0.015. (51) 

The range given is due to a variation of the top quark mass between 110 and 
170 GeV; however, we have not taken into account the uncertainties in the 
J'b determination from the lattice (±40 MeV) and the experimental error in 

X. 

Finally, we want to stress that the major part of our calculation as given 
in section 2 does not rely on vacuum insertion assumption. Eventually, 
the matrix elements of the four fermion operators should be evaluated non- 
perturbatively, e. g. by lattice calculations. Since lattice calculations work 
best at low scales, our formulae are needed to connect the low scales with 
the b mass scale where the matching to the high energy theory is done. 
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Figure 2: Diagrams with insertions of the local operators 
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